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A substantial magnetodielectric effect is often an indication of coupled magnetic and elastic order, 
such as is found in the multiferroics. However, it has recently been shown that magnetism is not 
necessary to produce either a magnetoresistance or a magnetocapacitance when the material is inho- 
mogeneous [M. M. Parish and P. B. Littlewood, Phys. Rev. Lett. 101, 166602 (2008)]. Here, we will 
investigate the characteristic magnetic-field-dependent dielectric response of such an inhomogeneous 
system using exact calculations and numerical simulations of conductor-dielectric composites. In 
particular, we will show that even simple conductor-dielectric layers exhibit a magnetocapacitance, 
and thus random bulk inhomogeneities are not a requirement for this effect. Indeed, this work es- 
sentially provides a natural generalisation of the Maxwell- Wagner effect to finite magnetic field. We 
will also discuss how this phenomenon has already been observed experimentally in some materials. 



INTRODUCTION 



Macroscopic inhomogeneities can have a surprising im- 
pact on the magnetotransport of a material. A spectac- 
ular example of this is the anomalous transverse magne- 
toresistance observed in doped silver chalcogenides [J . 
Here, the presence of inhomogeneities distorts the cur- 
rent paths in such a way that the Hall part is mixed into 
the longitudinal part of the response, thus giving rise to 
a large, linear magnetoresistance that persists across a 
range of magnetic fields and temperatures [H, 0| • Since 
this is a classical effect that does not require magnetism, 
it is relatively insensitive to temperature and thus poten- 
tially useful for technological applications. Indeed, inho- 
mogeneities have already been exploited in the design 
of magnetic sensors: semiconductor devices with "ex- 
traordinary magnetoresistance" rely on the presence of 
a metallic inclusion to distort current and increase the 
device's resistance in a magnetic field @ . 

There are similar non-trivial effects in the case of 
the ac dielectric response: Catalan has shown that a 
material can display a substantial magnetocapacitance 
when an intrinsic magnetoresistance is combined with 
Maxwell- Wagner extrinsic effects such as contact effects 
and bulk inhomogeneities [7[. Subsequently, Littlewood 
and I showed that a magnetocapacitance can be induced 
by inhomogeneities alone and magnetism is not even nec- 
essary Q. In particular, we found that a finite mag- 
netic field induced a characteristic dielectric resonance, 
though no inductive element was present in the system. 
As we discuss below, this is a counter-intuitive phe- 
nomenon that every experimentalist should be aware of 
when probing the dielectric response of a material, partic- 
ularly given that the existence of a magnetocapacitance 
is often taken as an indication of coupled ferroelectric 
and magnetic order. It is already known that compos- 
ites consisting of a ferroelectric and a ferromagnet can 
generate a magnetoelectric coupling in the absence of 



any intrinsic coupling [9]. Here, we examine how sim- 
ple conductor-dielectric composites can yield a dielectric 
response that depends on magnetic field in the absence 
of any magnetic or ferroelectric order. This effect is in- 
sensitive to the microscopies and is thus generic, being 
only dependent on the distribution of local capacitive and 
conductive regions. The only relevant material parame- 
ters are the static dielectric constant e of the capacitive 
regions, and the resistivity tensor p (or conductivity ten- 
sor a = p^ 1 ) in the conducting regions. As such, this 
magnetodielectric effect can occur in a variety of inho- 
mogeneous materials: for instance, it has already been 
observed experimentally in nanoporous silicon [lCt Jill 
and potentially in the manganite La2/3Cai/3Mn03 [12j . 
as discussed in Sec. [Vj 

In the following, we will investigate the characteris- 
tic magnctic-field-depcndcnt dielectric response of two- 
dimensional (2D) conductor-dielectric composites using 
exact calculations, numerical simulations, and the ef- 
fective medium approximation. In particular, we will 
show that even simple conductor-dielectric layers exhibit 
a magnetocapacitance, and thus random inhomogeneities 
in the bulk are not a requirement for this effect. Rather, 
it is sufficient to have conductor-dielectric interfaces that 
are perpendicular to the flow of current (or the move- 
ment of charges) . This also appears to be consistent with 
experiment, since there is an observed link between mag- 
netocapacitance and interfaces with free charges [l3| . We 
will also discuss how an intrinsic magnetoresistance could 
enhance this effect. 



II. MODEL AND METHODS 

We consider a classical 2D composite medium consist- 
ing of purely dielectric regions (defined by dielectric con- 
stant e) and purely conducting regions. The latter re- 
gions have the following resistivity tensor in a transverse 
magnetic field H = Hz: 
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where we have assumed that the conductor is isotropic. 
We focus on the simplest case where p xx — p and p xy /p — 
pH = j3, with p the carrier mobility. Locally, the current 
j(w,r) is related to the electric field E(w,r) via Ohm's 
law E(w,r) = p(w,r)j(w,r) = (iivi(ui, r)) _1 j(w, r), and 
globally the system is driven by electric field (E(w)) at 
frequency uj, where (...) corresponds to a volume average. 
The measured response averaged over the whole system 
is then the effective resistivity defined from (E(w)) = 
p e (w)(j(a>)). For the standard experimental setup where 
the boundary conditions are such that (j y ) = as 
in Fig. [TJa), the components of the dielectric function 
that are actually probed are then given by e xx {uj) = 
(iwpe^xM)" 1 an d e xy {iJ) = (iupe^uj))^ 1 . The lon- 
gitudinal response s xx (uj) is the usual dielectric response 
measured in experiment, while the transverse response 
e xy {oj) can be extracted from a measurement of the trans- 
verse electric field E y (uj) = E x (uj)s xx {lu) / e xy (uj). 

For the isotropic two-component media considered in 
Sec. IIV1 one can make use of the self-consistent effective 
medium approximation outlined in Refs. [l4T - fl6l |. Here, 
we first imagine that we have a single inclusion embedded 
in an effective medium and then we average over all inclu- 
sions to self-consistently determine the response of this 
effective medium. This amounts to solving the coupled 
equations: 
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FIG. 1: The basics of probing and modelling dielectric re- 
sponse. Diagram (a) depicts the general measuring setup, 
where a rectangular sample is subjected to an oscillating elec- 
tric field E x x and a static transverse magnetic field H. The 
response of a general 2D conductor-dielectric composite can 
be simulated with a network of four-terminal elements (b). 
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FIG. 2: Different configurations of conductor-dielectric lay- 
ered composites, where (j y ) = on the upper and lower 
boundaries, and we have periodic boundary conditions for the 
left and right boundaries. The interfaces between conductor 
p and dielectric e can be perpendicular (a) or parallel (b) to 
the current flow. 



where a e = p~ x and pi is the volume fraction of the i-th 
component. For the special case where the fractions are 
equal (pi = P2 — 1/2), we can exploit a symmetry trans- 
formation for the electric field and current density fll\ - 
[l9| to derive an exact result for the effective dielectric 
function of the composite medium Q. Note that the ef- 
fective medium approximation recovers this exact result 
for equal fractions, and becomes formally exact in the 
limits pi — > 0, P2 — > 0. Thus, we expect it to provide 
an accurate approach for investigating isotropic media 
across all volume fractions. 

For composites with anisotropic inhomogeneities like 
the conductor-dielectric layers in Sec. IIII1 one requires 
a more brute force numerical approach. In this case, 
we discretise the system into four-terminal elements as 
shown in Fig. [ljb) , where the voltages Vi at the terminals 
of each element are linearly related to the currents Ii 
via an impedance matrix: Vi = Zijlj. Note that each 
voltage is defined with respect to the element centre and 
incoming currents are defined as positive. In a purely 



conducting region, Z^- 



\(5ij+PMij), where 
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In a purely dielectric region, we simply have Z\j = 
(2iu]e)~ 1 6ij. To simulate a given layed composite such 



as in Fig. [21 we construct anfxM rectangular network 
of these four-terminal elements and then take the limit 
of large N and M (keeping N/M and the proportions 
of each phase fixed). In general, the response converges 
rapidly to that of an infinite network. 



III. CONDUCTOR-DIELECTRIC INTERFACES 

In this section, we consider simple conductor-dielectric 
layers (Fig. [2]), which constitute the simplest realisation 
of anisotropic inhomogeneities. They also allow us to in- 
vestigate the behaviour of conductor-dielectric interfaces 
and they can thus potentially describe the effect of con- 
tacts in experiment. Indeed, a magnetocapacitance has 
already been observed in non-magnetic Schottky barri- 
ers goj. 

In general, we find that the behaviour is strongly de- 
pendent on geometry. When the conductor-dielectric in- 
terfaces are parallel to the current flow as in Fig. [2f b) , 
the dielectric response does not depend on magnetic field 
and we simply obtain e xx (uj) = e + (iujp)^ 1 . By con- 
trast, when the interfaces are perpendicular to the flow of 
charge as in Fig.[2ja), we have a substantial magnetodi- 
electric effect. Firstly, at zero magnetic field, we recover 
the dielectric relaxation 2l[ expected from the Maxwell- 
Wagner effect — for a square unit with equal proportions 
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of each phase, i.e. L\ = 2Li 



L/2 in Fig.[2ja), this gives 

1 — ILOT 
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where the time constant r = pe. Here, at the character- 
istic frequency lot = 1, there is a rapid change in S[£ 1X ] 
and an associated peak in 3[£ xa; ]. Then, in the presence 
of a magnetic field, the peak shifts to smaller lot with 
increasing as shown in Fig. El Indeed, at large fields 
^ > 1, the characteristic frequency becomes Plot = 1 
and the behaviour evolves into a dielectric resonance, 
where the capacitance ^[exs] can become negative. This 
is similar to the magnetocapacitance of isotropic com- 
posite media [8[ considered in Sec. [IVJ We also obtain a 
similar magnctodielectric effect for other configurations 
of Fig. [Ha), but the shift of the peak position in ^[e^a;] 
often only approximately obeys Plot = 1. 

To gain insight into this effect, we consider a modi- 
fied version of the setup in Fig. [2ja) where we insert a 
perfectly conducting metallic interface between the con- 
ductor and dielectric. Now, if we ignore the dielectric 
for the moment and consider dc transport through the 
conductor, then such boundary conditions will give rise 
to a magnetoresistance 22]. In particular, if the conduc- 
tor is square {L\ = L), then the effective resistivity is 
exactly p*(H) = Py /l + P 2 0,111]. Thus, the response 
of the conductor-dielectric composite becomes (assuming 
L 2 = L/2): 
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In the limit /3 ^> 1, we essentially obtain Maxwell- Wagner 
relaxation with lot replaced by Plot. Thus, we see that 
the effect of the interface is to mix the Hall component 
into the response (both real and imaginary parts) and 
replace p with the Hall resistivity p/3. Note that we do 
not obtain a dielectric resonance here and thus it appears 
that this feature requires non-zero Hall fields E y within 
the dielectric itself, which is not the case for perfectly 
conducting interfaces. 



IV. ISOTROPIC COMPOSITE MEDIA 

We now turn to the dielectric response of a 2D isotropic 
two-component medium where we once again have purely 
capacitive and purely resistive regions. When there are 
equal proportions of each component, we have the fol- 
lowing exact results for the dielectric function Q 
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Note that e xx ~ l/v^ m both low and high frequency 
limits, in contrast to the layered case with equal propor- 
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FIG. 3: Magnetodielectric effect in conductor-dielectric layers 
for the geometry in Fig.^a) with Li = 2L 2 = L/2. For /3 = 
we have ordinary Maxwell- Wagner dielectric relaxation, while 
for /3 > 1 it becomes a dielectric resonance with characteristic 
frequency Plot = 1. 



tions. However, for /3 > 1, we have a dielectric reso- 
nance phenomenon similar to Fig. [3l with characteristic 
frequency Plot = 1. 

We can also allow for the possibility that the resistive 
component has an intrinsic magnetoresistance p X x{H) by 
considering the more general expression for the resisitiv- 
ity tensor in Eq. ([1]). In this case, we obtain 
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Here, we find that the size of the magnetocapacitance 
(i.e. the change of $t[e xx (Lo)]/e with P) depends crucially 
on the ratio p xy /p xx . Thus, we can have an enhanced 
magnetocapacitance for a negative intrinsic magnetore- 
sistance, where p xy /p xx ~ pP/p X x is large at finite mag- 
netic field. However, without a Hall component (i.e. 
Pxy = 0), we simply have s xx (lo) = e/y/iujep xx and we 
only have a magnetocapacitance when there is an intrin- 
sic magnetoresistance p xx {H). In this case, we have no 
relaxation or resonance phenomena. 



A. Effective medium approximation 

To investigate all volume fractions of the components, 
we use the self-consistent effective medium approxima- 
tion. Following Ref. [l6j | , we simplify the problem by first 
transforming to a frame where the resistivity tensors of 
each phase are scalar. Thus we take 



E' = E + bRpJ 
J' = cJ + Rp^B 



(9) 
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where R is a rotation by 90 degrees: 

R = ( i o 1 

and the transformation constants are 



(11) 
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using Eq. ([1} for the first component pi. Note that 
these coefficients are real, even though the dielectric 
phase (component 2) has a purely imaginary resistivity, 
P2 = (iujs)^ 1 . In the transformed frame, the (dimension- 
less) resistivities of the two components are 
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From Eq. the effective resistivity of the transformed 
composite is then 
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where p is the volume fraction of the dielectric phase. 
Transforming back to the original frame then yields the 
final result for the longitudinal effective resistivity: 
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Note that taking equal proportions (p = 1/2) yields the 
exact result described earlier — we can recover Eq. 
using e xx (oj) — {iLop e ^ xx )~ x . In what follows, we assume 
that we simply have p xx = p and p xy — ftp. 

For all volume fractions p, one can show that there 
is always a dielectric resonance at ftur = 1 for large 
magnetic fields Q. Moreover, we find that the dissipation 
at the peak only depends on p to leading order in ft: 



Pe 



l-2p+ v /l + 4p(l-p) 



■iO(l/ft) (14) 



where we see that p e , xx / p — > 00 when p — > 1, as expected. 
Contrast this with the layered case (e.g. Eq. ([5])), where 
we find that p e , xx /p at the peak increases with increasing 
magnetic field ft. 
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FIG. 4: Magnetocapacitance of isotropic conductor-dielectric 
composites at fixed frequency uit = 1 for different dielectric 
volume fractions p. In the limit of high magnetic field ft — > 00, 
&[e M ]/e ->■ and U\e xx ]/e -t 2p- 1 for p < 1/2 and p > 1/2, 
respectively. 



If we fix the frequency oj, we also find a sizeable mag- 
netocapacitance that depends on p. Figure @] depicts the 
magnetocapacitance for lot = 1, and we see that both 
the size and sign of ^t[e xx ]/e at large ft is determined by 
p. Indeed, in the limit ft 00, ^t\e xx \/e only remains 
finite when p > 1/2 and the dielectric phase percolates 
through the whole sample. 

Further insight can be obtained by deriving the limit- 
ing behaviour of the ac response from Eq. (|13p . Lets first 
consider the case where the magnetic field ft — 0. Then 
in the limit lot — > 0, we have 



p<l/2 
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Likewise, in the limit lot — > 00, we have 
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(16) 



Here, we can clearly see the effect of percolation, since 
the conducting regions dominate the transport below the 
percolation threshold of the dielectric (p < 1/2), while 
the dielectric regions dominate above (p > 1/2). In par- 
ticular, we see that the dependence on p is inverted for 
low and high frequencies, since the dielectric regions ef- 
fectively behave insulating in the former case and metal- 
lic in the latter (where there is insufficient time for the 
dielectric to build up charge). However, for the special 
case p = 1/2, i.e. at the percolation threshold, we simply 
have p e ,xx — pj \/iioT and thus the same behaviour in 
both limits. 

When there is a large magnetic field (ft ^> 1), we obtain 
the same result for lot —> 0, or ftooT <C 1, but there can be 
markedly different behaviour when ftcor > 1. Specifically, 
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cal transport in the semiconductor as being activated: 
T , where A is the activation gap. Then our 
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FIG. 5: Position of the resonance peak in the imaginary part 
of the dielectric response in terms of inverse temperature 1/T 
as a function of magnetic field H . Triangles are data for 
nanoporous silicon taken from Ref. while the line corre- 
sponds to a fit from our classical model of conductor-dielectric 
composites. 



for p — > oo and finite lot, we obtain 
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/3\Aljt 
1+iojT ' 
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p = l/2 
p> 1/2 
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which matches the large j3 behaviour in Fig. Note that 
for p < 1/2, the system at high frequencies behaves like 
one at low frequencies, Eq. (|15p . This implies that the 
dielectric region behaves insulating rather than metallic 
as one might normally expect at high frequencies. This 
is because at large magnetic field, the current paths are 
distorted such that they avoid any metallic region (cf. 
Ref. [a]) and thus any metal effectively behaves like an 
insulator. 



V. DISCUSSION 

The dielectric response investigated in this paper only 
depends on the presence of macroscopic inhomogeneities 
and can thus potentially be observed in a variety of ma- 
terials. There is already evidence to suggest that it 
accounts for the magnetodielectric effect in graphene- 
polyvinyl alcohol nanocomposite films [24| and in nickel 
nanosheet-Na-4 mica composites [25| . Furthermore, our 
results are possibly relevant to rough interfaces in het- 
erostructures [26(. 

Experiments on nanoporous silicon [ToL [llj have re- 
vealed a magnetic-field-dependent dielectric resonance 
that is consistent with our model. In this case, the di- 
electric response was measured at fixed frequency U) — 
100 kHz as a function of temperature T — 25-40 K 
and magnetic field H = 0-32 T. However, varying the 
temperature is equivalent to varying ujt if we approx- 
imate e as temperature independent and the electri- 



predicted resonance occurs at the temperature ksT = 
— A/ log(/3wepo)j so that the dielectric resonance shifts 
to higher T with increasing H , as indeed was observed in 
experiment. Assuming that /3 = 1 corresponds to H ~ 
IT, we can fit the data with uj = l/po£ - 4x 10 10 Hz 
and A ~ 30 meV, as shown in Fig. [SJ Note that the size 
of A is consistent with activation from an impurity band. 

Our effect may also be relevant to the dielectric reso- 
nance observed in the manganite La 2 /3Ca 1 / 3 Mn03 jus t 
above the ferromagnetic transition temperature [12j . 
However, here the situation is more complex since there is 
magnetism involved and the effective "composite" corre- 
sponds to the phase separation between magnetic metal 
and charge-ordered insulator. Moreover, the size of mag- 
netic metallic domains could change with magnetic field. 
This all requires further investigation: in principle, one 
might be able to probe the latter possibility by measuring 
the Hall component e xy , as discussed in Ref. Q. 

Finally, with regards to experiment, there is the ques- 
tion of when our classical model breaks down as the size 
of the sample is reduced. Formally, our approach is only 
valid when \jui is greater than any microscopic timescale, 
e.g., the scattering time of charge carriers within the con- 
ducting regions, and when the size of the inhomogeneities 
is greater than any microscopic lengthscale, e.g., mean 
free path. However, one may be able to extend our re- 
sults to smaller sizes using an appropriately redefined e 
and p. 



VI. CONCLUSION 

In this paper, we have investigated 2D conductor- 
dielectric composites and we have shown how these can 
produce a magnetic-field-dependent dielectric response 
without any magnetism. Using exact results and the 
effective medium approximation, we have derived ex- 
pressions for the dielectric response of isotropic compos- 
ite media and we have explicitly revealed the behaviour 
at the resonance (3ujt = 1 and in the limits ujt — > 0, 
ujt — > oo and /3 — > oo. Moreover, we find that an intrin- 
sic magnetoresistance could enhance this magnetodielec- 
tric effect from inhomogeneities. We have also performed 
numerical simulations of layered composites to show that 
conductor-dielectric interfaces play a key role in this phe- 
nomenon. Finally, we have discussed how our predicted 
magnetodielectric effect has already been observed ex- 
perimentally in some materials. An open question is how 
our results generalise to higher dimensions and this will 
be the subject of future work. 
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